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We study the self forces acting on static scalar and electric test charges in the spacetime of a 
Schwarzschild black hole. The analysis is based on a direct, local calculation of the self forces via 
mode decomposition, and on two independent regularization procedures; A spatially-extended par- 
ticle model method, and on a mode-sum regularization prescription. In all cases we find excellent 
agreement with the known exact results. 
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0\ • I. INTRODUCTION AND OVERVIEW 

^ I The problem of calculating the gravitational wave forms generated by compact objects orbiting black holes is of 
) ^ ' crucial importance for the detection and the interpretation of observations by gravitational wave observatories such 
as LISA [0. A major step towards the calculation of the wave forms is the computation of the gravitational radiation 

■ reaction forces, acting on the compact object. The generation of very accurate templates for the waveforms detected 
' from a system of a compact object in orbit around a supermassive black hole is an extremely hard task. It is likely 

that one would need to have accurate templates for as many as 10^ orbits. For such a system, accurate templates 
^ [ are necessary for detection, because the predicted signal to noise ratio for LISA is approximately of order 10 for a 1 
year integration time. Lack of accurate templates would results in a loss of a factor of roughly the square root of the 
number of orbits in sensitivity , which would result in signal to noise ratio below the detectability threshold. 

Several methods have been suggested for the calculation of radiation reaction. One approach follows Dirac's method 
for obtaining the Abraham-Lorentz-Dirac equation for an electric charge in arbitrary motion in Minkowski spacetime 
' 1^. In that approach one imposes local conservation laws on a tube surrounding the world line of the particle, 
I and integrates the conservation law across the tube, thus obtaining the equations of motion, including the radiation 
' reaction effects. In Dirac's approach the infinities which are related to the divergence of the particle's field on its world 
^H, line are removed by a simple mass renormalization. This method was used by DeWitt and Brehme [Q to generalize 
I ■ Dirac's analysis for a general curved background. More recently, Mino, Sasaki, and Tanaka used a similar method for 
^j^' the case of a massive particle coupled to linearized gravity Q . Quinn and Wald Q formulated recently an axiomatic 

■ approach for the calculation of radiation reaction. In that approach, the infinities are removed by comparing the forces 
. !^ I in different spacetimes. However, at present it is unclear how to apply the Quinn- Wald formal approach directly to 

■ practical calculations. The main difficulty arises from the calculation of the "tail term", which is difficult to compute 
Jh [ even in the slow motion, weak field limit |@,||. 

Another approach is based on arguments relating to the balance of quantities, which are constants of motion in 
the absence of radiation reaction, specifically the energy and angular momentum. Such balance arguments involve 
integration of the flux of an otherwise conserved quantity over a boundary which consists of a distant sphere and 
the horizon of the black hole |]9|-|l6|]. Although these methods are quite successful to very high relativistic order 
in Schwarzschild spacetime, they are problematic for the more interesting problem of motion in the spacetime of a 
spinning black hole, because of an inherent difhculty: For the problem of motion in the Schwarzschild spacetime the 
motion is completely determined by the rate of change of the energy and the angular momentum, which are additive 
constants of motion in the absence of radiation reaction. However, for general orbits in the spacetime of a Kerr 
black hole there is a third constant of motion, i.e., the Carter constant. The Carter constant is non-additive, and 
consequently it cannot be obtained by methods which are based on balance arguments. (Such methods can be used 
for circular and equatorial orbits around a Kerr black hole, because then the evolution of the Carter constant is trivial: 
it is given completely by the evolution of the energy and the azimuthal component of the angular momentum.) In 



addition, such methods suffer also from other difficulties |17|: they usually yield only the time average of the radiation 



reaction force, such that for any fast evolving system they would be inherently inaccurate. In addition, they fail at 
obtaining the conservative piece of the radiation reaction force. 

A different approach for the calculation of the gravitational radiation reaction is based on a direct, local calculation 
of the self forces acting on the compact object. Obviously, knowledge of the instantaneous forces acting on the orbiting 
object would allow for the calculation of the orbital evolution. Such a direct approach for the calculation of the self 
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force was suggested by Gal'tsov ||T^. However, Gal'tsov's approach is based on the radiative Green's function (i.e., 
the "half-retarded minus half-advanced" potential), which is not causal in curved spacetime, because it requires the 
knowledge of the complete future history of the object in motion [|9|. A causal approach, which is based on the 
retarded Green's function rather than on the radiative one and consequently is more in the spirit of relativity theory, 
is much more desirable. 

Recently, a local approach for the calculation of the self force, which is based on the retarded field, and on a 
Fourier-harmonic mode decomposition of the field and the self force, has been proposed |p^-pl[. This approach has 
two very important advantages: First, when the field is decomposed into modes, each mode satisfies an ordinary 
differential equation rather than a partial one, and consequently the solution for each mode is considerably simpler. 
Second, and most importantly, each mode of the self force turns out to be finite. Indeed, the total self force, which 
is obtained when one sums over all modes, very frequently diverges, but this difficulty is met only at the summation 
over all modes step: The treatment of the individual modes is free from divergences. 

This approach was used by Ori ||l9|,0 for the calculation of self forces acting on objects in motion around a black 
hole. Ori suggested a regularization prescription which is based on the assumption that the divergent piece of the self 
force is proportional to the four-acceleration of the charge. One can then use a simple mass-renormalization procedure 
by re-defining the mass of the particle to include the divergent piece. The self forces in general is expected to diverge. 
Therefore, a crucial ingredient for the calculation of the self force is the regularization method which one uses. We 
note that in any regularization prescription in the gravitational case (i.e., when the particle has a non-zero mass) one 
faces a gauge problem. This difficulty, however, does not arise in the cases of scalar or electromagnetic fields, because 
the force in these cases is gauge independent. Therefore, consideration of scalar or electromagnetic charges is of some 
value, as they correspond to easier cases, where much of the difficulties related to self forces are already present, yet 
one does not have to solve also the gauge problem. 

The self force is expected to be proportional to the particle's (charge)^. The correction to the orbit is therefore 
also of order (charge)^, and the correction to the self force is consequently of order (charge)"*. When the charge of the 
particle is much smaller than the mass of the black hole, this correction is negligible. In this paper we shall study the 
self force only to leading order, i.e., we shall study the self force to order (charge)^. 

We present in this paper two independent regularization procedures for the self force, which are successful for the 
problem of static charges (both scalar and electric) in the spacetime of a Schwarzschild black hole. (These procedures 
were also found to be successful for the regularization of the radial component of the self force for scalar or electric 
charge in uniform circular motion in flat spacetime pT| , p2[ ). We hope that similar methods (or their generalizations) 
would be relevant also for more complicated and realistic problems, e.g., the self forces acting on a compact object in 
circular motion around a Schwarzschild black hole, and ultimately, the self force on a compact object in motion in a 
generic orbit around a Kerr black hole. 

Let us consider first spatially-extended particles (we still assume that the extension of the particles is smaller that 
the typical radius of curvature and the typical scale of inhomogeneity of the field). The divergent piece of the self 
force, in addition to being proportional to the four-acceleration, is also expected to be inversely proportional to the 
spatial extension of the particle. In the limit of a point-like particle, this is the source for the divergence of the 
force. One should therefore be able to obtain a regularization procedure by considering a spatially-extended model 
for the particle, and then consider a sequence of smaller and smaller particles. The force acting on the particles would 
increase like the inverse of their size, and by removing this piece of the force one can expect to obtain the regularized 
self force, which is independent of the assumed internal structure, in the limit of vanishing spatial extension. The 
question of whether the regularized force depends of the way the point-like limit is taken is still an open question. 
A similar approach was used by Ori, who calculated the self forces acting on static scalar and electric charges in 
Schwarzschild and on the axis of Kerr black holes |2^. In Schwarzschild, Ori used a dumbbell model, where the axis 
was aligned either radially or tangentially. In Kerr, Ori used a radially-aligned dumbbell model. Whereas we use a 
mode-decomposition approach, which does not depend on the availability of an exact solution, Ori used the exact 
solutions for the scalar field or the electric potential, which are available for these cases, in order to calculate the self 
forces. 

We also consider a second, independent, regularization prescription. We consider a point-like particle. In that 
case the sum over modes is expected in general to diverge. Ori has recently suggested a mode-sum regularization 
prescription (MSRP) for the self force j2^j2^. Although MSRP is not fully developed as yet, it has already been 
shown to be valid for simple cases, such as a static scalar charge outside a Schwarzschild black hole. MSRP can 
possibly be generalized also for more complicated cases, such as massive particles in orbit around a Kerr black hole. 
If robust, MSRP can be of great importance for the generation of templates for the detection of gravitational waves 
from compact objects in motion around supermassive black holes. 

The organization of this paper is the following. In Appendix A we describe very briefly the main ideas of MSRP, 
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applied for a scalar charge in Schwarzschild. In Section || we discuss the self force acting on a static scalar charge in 
Schwarzschild spacetime. The result has been obtained by independent methods: For a minimally-coupled massless 
scalar field the self force vanishes |2^Jl^,^ . It is our approach which is novel: Our calculation is based on a direct 
computation of the self force mode by mode, followed by a summation over all modes, and finally on two independent 
regularization procedures. One regularization procedure is based on a spatially-extended particle model. We then 
consider the forces acting on a sequence of such particles with decreasing spatial extensions, and remove the divergent 
piece of the self force by a simple mass renormalization procedure. The other regularization procedure is based on 
MSRP. We find that both methods are successful in obtaining the correct result. In Section |l| we consider the 
analogous problem of the self force acting on a static electric charge in Schwarzschild spacetime. Also in this case, 
the result is not new. This problem has been considered by several authors: DeWitt and DeWitt calculated the 
radiation damping forces (both nonconservative and conservative) acting on a slowly moving electric charge in the 
far-field regime, and found that there was a repelling self force, which lowered the much stronger gravitational pull of 
the black hole, and made a retrograde contribution to the periastron precession. Vilenkin considered the electric 
charge to be very far from the black hole (specifically, he assumed the position of the charge to be at rp ^ M, where 
M is the mass of the Schwarzschild black hole), and again found that there was a repelling conservative self force. 
Smith and Will |2|] and Frolov and Zel'nikov [ ^0|j2^ ] were able to solve for the force exactly, for all positions of a 
static charge in Schwarzschild spacetime, and found that the repulsive radial self force was f~^^^^ — e^M /r^ (in the 
frame of a freely falling observer who is instantaneously at rest at the position of the charge). Also in this case of 
a static electric charge we present a direct approach for the calculation of the self force, which is based on mode 
decomposition, summation over all modes, and force regularization procedures similar to those we apply in the scalar 
case. In Section IV we summarize our methods and results. 



II. STATIC SCALAR CHARGE 
A. Mode decomposition of the force 

Consider a point-like scalar test charge in the Schwarzschild spacetime, held fixed by some external force. Our 
aim here is to calculate the contribution of the self force to the total force needed for keeping it fixed. The result is 
well known |26 1^J2^: the contribution of the self force to the total force vanishes. The linearized field equation of a 



minimally-coupled, massless scalar field $ in the Schwarzschild geometry, which is described by the line element 



ds^ +[l dr^ dif. 



dt^+il 

r J \ r 

where dQ^ — dO^ + sin^ 6 d^p^ , is given by 

V^V^$(a:") = -4^p(a;"), (1) 
where denotes covariant differentiation, and where the charge density 

P^.r^r^^l^^l^SM. (2) 



Here, q is the charge, r is its proper time, and g is the metric determinant. The mass of the black hole is denoted by 
M . The world line of the charge is given by xI^{t). In what follows we use the usual Schwarzschild coordinates: The 
radial Schwarzschild coordinate is defined such that spheres of radius r have surface area 47rr^, and t is the proper 
time of a static observer at infinity. We take the charge to be on the equatorial plane at the coordinates r = rg, (^s = 0, 
and 6 — 7r/2, without loss of generality. (Because of the symmetry of the Schwarzschild geometry the coordinates 
9 and tp can be rotated such that these would be the coordinates of any static charge at r — tq.) Because of the 
staticity, the scalar field is independent of the time, and we can decompose it into modes according to 



<I>(r,0,^)=£ ^ 0'(r)r'™(0,^), (3) 

/=0 m=-l 

such that the LHS of the wave equation (|^) is given by 
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The charge density is similarly decomposed into modes according to 



p = q- 



J{r - rg) 1 
3 u*{ro) 



(4) 



(5) 



;=o m=-l 



where it" is the four velocity of the charge, and a star denotes complex conjugation. We thus find the radial equation 
for 0'(r) to be 



2M\ ,, 2 , , , 1(1 + 1) ,, 



(6) 



To solve this equation we transform to dimensionless harmonic coordinates, i.e., we define f = (r — M)/M. In the 
harmonic gauge the radial equation is nothing but the Legendre equation^. We choose the two independent solutions 
of the corresponding homogeneous equation to be Pi{r) and Qi{r). The former is regular for 1 < r < f^, and the 
latter is regular for f > fo- (Note that the horizon of the black hole is located at f horizon = !)• The summation over 
all modes m is readily done, and we thus write the field at the point (r, 9, ip) due to a scalar charge q at the position 
{rs,0s,'Ps) as 



_<l_ 

M 

Pi 



9M 

1 ^(2? + l)Pi(cos7) 
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- M 



M 



e(r, - r) 



(7) 



Here, COS7 = cos^cos^g + sin0sin6'5 cos((/5 — ips)i and 0(x) is the Heaviside step function, i.e., Q{x) = 1 for a; > 
and Q{x) = for a; < 0. This solution for the scalar field $ is regular both at the black hole's event horizon and at 
infinity. In what follows as choose the angular coordinates such that both the origin and the evaluation point of the 
field lie on the equatorial plane, such that cos 7 = cos(v5 — tps)- 

The force which a scalar field ^E" exerts on a scalar charge q' is given by /„ — q' + w^w'''^',/?), where the four- 
velocity m'" is that of the charge q' . The scalar field ^ can be any scalar field, in particular the self field of the charge 
in question itself. In a sense, this is the scalar field analog of the electromagnetic Lorentz force [Q. Because of the 
staticity of our problem the only component of u'" which does not vanish is the temporal component. However, the 
temporal derivative of the field vanishes, and consequently the force is given only by fa = q'^ ,a- Because for scalar 
fields partial derivatives equal covariant derivatives, this is, in fact, the covariant equation for the force. Consider now 
two scalar charges, qi at (ri, 7r/2, ipi) and 52 at (r2, 7r/2, (^2)7 and r2 > ri. The force that qi exerts on (72 is given by 



^ That this should be the transformation is most easily seen from the following consideration. In the dimensionless coordinate 
X = r/(2M) the homogeneous equation is (1 — x)x(j)"(x) + (1 — 2x)4>'{x) + l{l + l)(j){x) — 0. This is a hypergeometric equation 
of the canonical form x{l — x)<j)" + [c — [a + b + l)x](j)' — abcf) — 0, for a = —I, b = I + 1, and c = 1. As 1 — c = c — a — 6, we 
know from the theory of hypergeometric functions that the homogeneous equation can be transformed to Legendre's equation. 
The variable of the hypergeometric equation x is then related to the variable of the Legendre equation by the transformation 
a; = (1 + f)/2. In view of the definition of x, we find that f is nothing but the dimensionless harmonic coordinate. We are then 
assured that transformation to f would yield Legendre's equation with solutions Pi(f) and Qi(f) |p4|. 
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Similarly, the force that q2 exerts on qi is given by 



A'P V f2 



/=0 



r , M ,f'ri~M\ fr2-M\ , , 

X P,[cos(^2 - Vi)\P[ ( ) \j^jr~ ) 



Here, a prime denotes differentiation with respect to the argument. Each of the force components is evaluated at 
the position of the charge on which the force is exerted. We are interested to calculate the self force acting on a 
point particle. Namely, we are interested in identifying q\ with qi- When this is done, one naturally finds that the 
total force diverges (although each of the / modes of the force is still finite). Next, we describe two regularization 
procedures for the self force acting on a point particle, which yield the desired result. 




B. Regularization procedures 

1. Extended particle model: Inclined dumbbell 

A well-known classical renormalization scheme is to consider a spatially-extended particle model, and then consider 
the limit of vanishing spatial extension, in the spirit of the classical Abraham-Lorentz-Poincare electron models. 
However, as is well known point- like particles are problematic in General Relativity even to a greater extent than 
they are in electromagnetic theory because of the non-linearity of the Einstein equations |3^ ]. Still, in some sense, 
one can be hopeful that as the particle becomes smaller and smaller, the deviation of its world line from a geodesic 
becomes insensitive to the particle's internal structure. The simplest particle model is a dumbbell model, consisting of 
two point-like charges at the two edges of an uncharged rigid rod, whose length is smaller than the typical scales of the 
inhomogeneities of the gravitational and scalar or electric fields (having in mind that we shall later consider the limit 
of vanishing spatial extension). Although this is a very simplified model for a particle, it can be simply generalized 
to more realistic models, bearing in mind that a general extended (classical) object can be construed as comprised 
of many point-like particles, and the self interaction of a general extended object can be obtained by summing the 
contributions of all pairs of point-like particles, each pair being, in fact, a dumbbell. We shall thus model the particle 
to be a dumbbell, with two equal charges qi = q2 = e/2, e being the total charge of the particle. Because of the 
symmetry of the geometry, the simplest configuration is to align the dumbbell axis in the radial direction. That way, 
we still maintain axial symmetry, and the dumbbell axis is aligned along a geodesic. However, we shall see below that 
despite the fact that with a radial dumbbell axis one indeed recovers the known and correct result for the self force, an 
important feature of a general extended particle model is missing, specifically, the mass-renormalization aspect of the 
force regularization procedure. This happens because the coefficient of the divergent piece of the bare force vanishes 
if the alignment of the dumbbell is radial (in the scalar case). As we are interested primarily in the regularization 
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procedure, we shall consider here a more complicated case, where the dumbbell is not aligned radially. Consequently, 
we shall take the dumbbell axis to be inclined at some angle from the radial directionR. 




FIG. 1. The geometry of the charge splitting in the equatorial plane: The charge e is split into two half charges, gi = g2 = e/2. 
The charge q\ is placed at ro — e and the charge q2 is placed at ro + e. The total angular separation between q\ and q2 is 2A<^. 
The horizon of the black hole is at r = 2M. 

Specifically, we take r2 — tq + e, ri — tq — e, ip2 — Aip, and ipi = —Aip. For concreteness, we take Aip = ae, 
such that when we make e smaller, we also reduce Atp proportionally, and we take e <C 2M. Figure ^ illustrates the 
geometry of the charge splitting in the equatorial plane (recall that the coordinates can always be rotated such that 
the splitting is in the equatorial plane) : The charges qi at ro — e and q2 at + e are also separated angularly by an 
angle of 2A(p. When wc take the limit e ^ we simultaneously take the limit Aip ~> too, such that the point-like 
charge is located at the intersection of the circle of radius ro and the bisector of the angle between qi and q2 ■ 

Let us consider only the radial force which acts on the dumbbell. (Because the acceleration is purely radial, we 
expect only the radial component of the self force to diverge.) The total (bare) self force which acts on the dumbbell 
is made of four contributions. Schematically, 

= + + + (12) 

f^^ being the radial component of the force which the charge qi exerts on the charge qj . Let us consider this force in 
the point-like particle limit. Now, Z*"'**' is the self force on a point-like scalar charge e. However, both f^^ and /^^ 
are equally well the self forces on point-like scalar charges, which are identical to the original charge e in all respects, 
except for the fact that they each have charge e/2. As the self force is proportional to the charge squared, it implies 
that = = /t°tai/4. Consequently, 

Because we need to sum vector components, we have to perform the summation at a common point, which for 



mathematical complication occurs if we take the dumbbell axis to be in the d/ dip direction. Specifically, in that case the 
series expansion for the scalar field indeed converges, but not absolutely. Consequently, one is not allowed to differentiate term 
by term to obtain the force. This difficulty can most easily be illustrated in flat spacetime, where it already occurs: The electric 

scalar potential due to a static unit point charge is given hy V = X^^q -TTTfi(cos7). (See Ref. for details.) When the 

splitting is tangential, r< = r> = r, and the potential is given simply hy V — (1/r) P;(cos7). It can be easily checked 

that this series converges (although very slowly), but because of the oscillations it does not converge absolutely. When the 
radial positions of the source and the evaluation point are not equal there is an additional attenuation, thanks to which the 
series converges absolutely. 
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symmetry we choose to be (rp, 7r/2, 0)|^ Specifically, we need to transport the forces f^{r2) and /^^(^i) to (ro, it/2, 0) 
parallelly. For non-zero inclination angles the two edges of the dumbbell are not separated by a geodesic of the 
background geometry. The final result for the self force should of course be independent of the artificial spatial 
extension we assume (i.e., independent of the internal structure of the particle), of the parallel-transport route, of 
the point where we sum the forces, and of the specific way at which we take the point-like limit. It is still an open 
question whether the final result depends on the way the limit is taken. One might be worried about the introduction 
of ambiguities due to the arbitrariness in the choice of the parallel-transport route. Any ambiguity is of the order of 
the area enclosed by the two routes we compare, times the curvature. The area is of order e^, and the curvature of 
order Al/r^, such that the ambiguity is of order e^M/r^ . Because the Coulomb components of the individual forces 
cancel (see below), the leading order term in the total force is of order e~^, such that the ambiguity in the total force 
is of order e, and vanishes in the limit e ^ 0. That is, the final result is independent of the parallel-transport route. 
We note that the fact that the two edges of the dumbbell are not separated by a geodesic is not a problem of principle, 
because for a general extended object all pairs of the object's atoms interact, and most of them are not separated by 
geodesies. 

We perform the parallel transport of /^^^(r2) to (rp, 7r/2, 0) in two steps: first, along the radial route (r2, 7r/2, (^2) 
(ro, 7r/2, (^2), and then along the tangential route (ro, 7r/2, (^2) (''0, """/S, 0). Similarly, we parallel transport P^iji) 
from (ri, 7r/2, (pi) first radially to (ro, 7r/2, t^i) and then tangentially to (ro,7r/2,0). Note, that although we are 
eventually interested only in the radial component of the self force, we need, in fact, to parallel transport both the 
radial and the tangential components of the forces in the first sections of both routes, because when one parallel- 
transports a tangential component of a vector tangentially, it acquires a radial component (already in flat space). 
Another point to be made concerning the parallel transport is the following: The individual forces can be expanded 
in a power series in e, where the leading term is proportional to e~^. The self force is of order unity, i.e., of order . 
Therefore, one needs to perform the parallel transport accurately at least to order e^. The parallel transports along 
the radial routes is done as follows: The change in a covariant component of a vector in parallel transport along the 
d/ dx^ direction satisfies 5Va = V^pV^dx^, where F^^ are the connection coefficients, which equal the Christoffel 
symbols of the second kind in General Relativity. For the Schwarzschild geometry we find 5fr = Tl.^{r)frdr and 

= rf^{r)f^dr. Consequently, S(logfr) = d (log 

sections of the parallel transports /;^™ = ff'^\l IzIm'/ ' ^ ^ ^^'^ /"™ - -"p 
parallel transportation, the routes are tangential, such that 5fr = ^'^r{''')fv '^'P- (We do not need to find the change 
in the tangential component of the force as we are interested eventually only in the radial force.) That is, we need to 
integrate dfr = {f<p/r) dip. This can be done straightforwardly, and we find 



l-2M/r 

7, and 



and (5(log/^) = d(logr), such that in the radial 



joid(-^ncw/^oid) ^j^g second sections of the 




Pi{cos2Aip) 



/ 1 / 2M 
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\ r2 V n 



X ^(2/ + l)Pi 



1=0 




X [Pi(cosA.^) - P,(cos2A(^)]} 



(14) 



In Appendix B we describe briefly the numerical method we use for the evaluation of the series. We evaluate this 
force for various values of e (recall that we take Aip to be proportional to e). We find that /*°*''' diverges like e~^, 



^Note that the dumbbell is not symmetric about this point, because the invariant distances from this point to the two edges 
are not equal. However, there is no particular need for a symmetric dumbbell, and therefore we choose the model which is the 
simplest mathematically. 
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for very small values of e. This is indeed the expected behavior for the bare force. Classical mass renormalization 
can be used for the regularization of the bare force. Specifically, the divergent piece of the force is expected to be 
proportional to the acceleration, such that it can be absorbed in the mass of the particle. We use the exact solution 
for the scalar field IITI 



$ = q^l - ^ [(r - M)2 - 2(r - M) (ro - M) cos 7 + (ro - Alf - sin^ 7^ 



(15) 



and sum the mutual forces of the two charges at the dumbbell's edges at a common point, in the same way as above. 
Then, we expand the total force in a power series in e, where the leading-order term is of order e~^. This leading-order 
term is given by 



Jr 



2... , 2M\ 1 r 2ri^(ro-M) 



1 



- <; a- 

-3/2 



Af3 



X a 
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+ 1_ 2M 



2 / 2M\ 



,,r\ 3/2 



1 _ 2M 
'"0 



(16) 



Notice that the mass renormalization term depends on the value of the parameter a. That is, for different choices of 
a we re-scale the mass by a different quantity. However, the renormalized, physical self force is independent of a, as 
should indeed be the case. We note that this mass renormalization procedure does not depend on the availability of 



an exact solution. Below, in Section III, when we discuss a similar mass renormalization for an electric charge, we do 
not use the exact solution (although it is available). Instead, we use Eq. (^7|) for calculating the divergent piece of 
the self force. Even in cases where an equation analogous to Eq. ( |3^ ) is not available, the regularization procedure 
can still be done. In such a case one can extract the asymptotic divergence of the force at small separation distances 
from the bare force, by finding the asymptotic growth rate, and remove this piece from the bare force. 
We define the renormalized self force to be 



rcn /"total 



/■rcn _ £ 
J r — J 1 



f: 



div 



(17) 



The value of /™" is of course a function of e, and we need to take the self force in the limit of e — > 0. We find that 
the larger a, the greater is the number of modes over which we need to sum until f^™ converges and the oscillations 
are damped. Figure ^ displays the behavior of the sums over modes up to a certain value of the mode number I 
as functions of I for several values of the inclination parameter a. It is clear from Fig. |2| that for large inclination 
parameters one needs to sum over many modes. In addition, we also find that, with fixed a, the number of modes one 
needs to sum over scales like e"""^. When these two effects are combined, one finds that it is very costly numerically 
to consider nearly tangential splittings. 
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FIG. 2. The behavior of the sum over modes of the renormahzed force as a function of /, for different values of the inchnation 
parameter a. For all cases we take ro — 2.1M, and Aiy3 — 0.1. Top panel (A): a — lOAf ~^ (corresponding to e = 1 x 10~^i\f). 
Middle panel (B): a = lO^M"^ (corresponding to e = 1 x 10"^M). Bottom panel (C): a = lO^M"^ (corresponding to 
e = 1 X lO^^'Af). 

Figure || shows the renormalized force, i.e., = /*°*^' — /^'^ as a function of e for a non-zero value of the 
inclination parameter a. Similar results were obtained also for other values of a (but the numer of modes we needed 
to sum over depended, of course, on the value of a). The figure shows that for small spatial extension (small values of 
e) the renormalized force is linear in e, such that in the limit of vanishing spatial extension the force would equal zero. 
Notice that we can see deviations from the linear law for large spatial extensions. These deviations are expected, 
because the renormalized force, when expanded in a power series in e, contains contributions from all non-negative 
powers of e. The self force is the force on a point-like particle, i.e., the force in the limit e 0. Consequently, for any 
non-zero value of e we have contributions also from all positive values of e, which are dominated by the linear term in 
e for small values of e. 



el M 

FIG. 3. The renormalized self force f^"^ as a function of the spatial extension e. The charge is located at ro = 2.1M, and 
we choose a — 0.1. We sum the / modes here up to / = 2.8 x 10^. 



In the special case where the alignment of the dumbbell axis is radial (a = 0) , we find that the divergent piece of 
the force /^''^ vanishes identically, such that /*°*'^' is already renormalized. In this case we can sum the series in /*°*^' 
analytically and find the self force exactly. In fact, for any non-zero e we find for a radial dumbbell axis 
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(18) 



where we used the summation X];^o(^^ + '^)[Pi{^)Q'i{v) + Pi{^)Qi{v)\ = for y > x > 1. 

We thus find that if one considers a spatially-extended particle model for the particle, one can obtain a finite result 
for the self force in the limit of vanishing spatial extension, after performing a simple mass-renormalization procedure, 
which agrees with the well-known exact result . 



2. Mode-sum regularization 

In this section we use MSRP in order to find the self force on a point-like static scalar charge in Schwarzschild. 
MSRP is described briefiy in Appendix A, where the notation and definitions of the MSRP parameters are given. We 
note that our discussion here serves a dual purpose: First, it applies MSRP for a specific case, and obtains non-trivial 
physical results. Seconds, because our results can be compared with the final results for the self forces, which are 
already known, it predicts values for the MSRP parameters, which can then be tested analytically. 

In the case of a point-like particle, we find from Eq. that the bare force is given by 

barcp _ 1 [~2M 

Obviously, when this series is naively summed, the bare force diverges. In order to check the applicability of MSRP 
we first observe numerically that the / modes of this forces fl approach a non-zero constant as Z — s- oo, which we 
denote by f!j!° . Figure ^ shows the convergence of the I mode of the bare force to a constant, as ^ — > oo. The 
top panel of Fig. ^ shows ^^"^/' as a function of Z, for the first few values of the latter, and the bottom panel shows 
the difference between two consecutive I modes of the force as a function of I. We find that this difference scales like 
for large values of I, which implies that indeed the series converges to a constant. This behavior implies that 
the MSRP parameters and Cr vanish. (A non-zero value of Cr implies that the difference between two consecutive 
modes should scale like l^"^ ■) 
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Because the / modes of the bare force approach a non-zero constant as / ^ oo, it is clear that the sum over all 
modes diverges to infinity. That is, the source for the divergence comes from the contributions of the large-/ modes. 
Let us assume now that this divergence can be regularized by removing the large-/ contributions. That is, we assume 
that the large-/ contributions to the regularized self force die off with /. The only sensible way to do that is to 
subtract the asymptotic value of the modes (as / oo) from all the modes of the bare force. Although this procedure 
yields a finite result for the self force, it is not a priori clear whether that is the correct, physical result. However, 
because the self force is already known, this can be checked, and we can predict a value for a possible finite additional 
term for the regularization procedure, which can be then be tested analytically using MSRP. For obvious practical 
reasons, we do the summation over the modes only up to a finite value of /. We denote the approximations of the bare 
and regularized forces (which are obtained by summing over a finite number of modes) by '^''''°F/(ro) and ™'^i^/(ro), 
respectively. Then, we represent ^'^''^ f^^ (ro) by the /' mode of the force, for /' much larger than the / up to which we 
sum the series. In practice, we find that /' « 3/ suffices to a very good accuracy. Specifically, 
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(20) 



Figure |^ shows the bare force '^"°F/(ro) and the renormalized force ™"F/(ro) as functions of /. The bare force of 
course diverges for large values of /. However, figure jH] implies that the renormalized force vanishes for large / like 
Recall that the self force for this case is already known to be zero [^,|l^. Consequently, we infer that the value 
of the possible additional term for the regularization procedure is zero. Indeed, MSRP yields for this particular case 
dr = 0, which agrees with our result (see Appendix A). Because dr = 0, the regularized self force is given by MSRP 
to be simply 



t-iF, = ^['^-°/;(ro)-&.(ro)] 

1=0 



(21) 



where br 



bare 



fr°°j the regularization procedure is reduced to subtracting the asymptotic value of the modes 



of the bare force from all its modes, and then summation over all the modes. 




FIG. 5. The bare force and the renormalized force as functions of the / up to which we sum over the modes. Top panel 
(A): ''="^'=F^(ro) as a function of I. Bottom panel (B): '■™i?,'.(ro) as a function of /. For the renormalization procedure we use 
/' = 4.5 X 10*. The scalar charge is located at ro — 2.1M. 



We can also check the prediction of MSRP for the exact value of br- Recall that in this case br = — [g^/(2r'^)](l — 
M/r)/(l - 2M/r), and that, with = 0, MSRP predicts '"'^''Vr ^ 6r as / ^ oo. Figure | displays the difference 
between fl and br as a function of /. This difference behaves like /~^ for large values of /. This asymptotic behavior 
again implies that ar — and Cr = 0, as we found above. For rg = 2.1M, we find this difference to be 1.39 x 10~^ for 
/ = 4 X 10''. This agreement between the analytical prediction for br and the value to which the modes of the bare 
force approach at large values of the mode number provides a strong support for the validity of MSRP. 
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FIG. 6. The quantity ^'^'^'^ fr{ro) — &i (ro) as a function of I. The scalar charge is located at ro — 2.1M. 



III. STATIC ELECTRIC CHARGE 
A. Mode decomposition of the force 

An interesting case to study with our method is the case of a static electric test charge in Schwarzschild spacetime. 
This is interesting because it is known that in this case the radial self force does not vanish. This can give us two 
benefits. First, we can see whether our method can yield a correct non-zero result (a zero result cannot reveal a wrong 
factor, say), and second, we can use the exact expression for the result to evaluate the error in our calculation. The 
exact result for the self force in this case was found by Smith and Will |^ and by Zcl'nikov and Frolov [Q. The 
field of a static electric charge in the Schwarzschild spacetime was found in terms of a series expansion solution by 
Cohen and Wald H (see also ^^.) 

The Maxwell equation in curved spacetime are given by 

V^F'"' = Anf (22) 

where the Maxwell field strength tensor is given in terms of the four- vector potential by F^^^ — A^,^ — and where 
= puf^ is the four-current density, p being the charge density. Because of the staticity of the problem (both the 
charge and the fixed background geometry are static), all spatial components of both the vector potential and the 
current density vanish. The temporal component of Eq. ( p2| ) becomes 

^. , = -47rj*. (23) 
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In Schwarzschild coordinates this equation is explicitly written as 
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We next assume a series expansion of the form 

At (r, 0, ^) = 5] ^ i?' (r)y '™ {e,ip) (25) 
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and decompose the current density into modes 
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This current density corresponds to a total charge q, as is evident from 
We thus find the radial equation to be 
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= 47rg<5(r-ro)r'™*(-,0). 



(28) 



The basic functions which solve the corresponding homogeneous equation, with a convenient choice of normalization, 
are given by |3^,pql 
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and R^{r) — 1. The Wronskian determinant of these two basic solutions is Wi{r) — —{21 + \)/r'^ . The solution 
of the inhomogeneous equation (|2|) is thus 



-^H^) = „r , N -R; (^)0(^o - J') 
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where 



/(ro) = 47rglr'™*(?'0). 
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The function Ri{r) is regular both at infinity and at the black hole's event horizon. The summation over all modes 
m is straightforward, and we find that the I mode satisfies 
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For the monopole term (/ = 0) we find A^ = —{q/r)Q{r — rp) — (9/»'o)0(ro — r). Also in this case an exact solution 
is known |l^ , which is 



At = ^ 
ror 



Af - 



(r - Af )(ro - M) - Af 2 cos 7 



(r - Af)2 - 2(r - Af)(ro - Af) cos 7 + (ro - Af)2 - M"^ sin^ 7 



(34) 



The total covariant temporal component of the four-vector potential is obtained by summing over all I modes. The 
expression we thus find for At is identical to the expression given in Ref. and Ref. 1 17 . For the calculation of the 
force we need only the gradient of At with respect to r, which we simplify with the differential equation which the 
Legendre functions satisfy. We find that 



13 



, 9nr ^9 ir-2M){ro-2M) 



E 



2/ + 1 



P, 



, fr p-M 
M 



+ PI 



Q;(^)e(r-ro) 



0: 



e(ro - r) 



Pi (cos 7) 



g ro - 2M 

M2 ; 



^(2Z + l)Pz(cos7) 



+ Pi 



M ) \ M 
r - M\ /ro - M 



M 



Q'l 



M 



e(ro - r) 



(35) 



We note that we did not include in this expression the terms proportional to a delta function for the following reason. 
When the field is evaluated at any point which is not the position of the charge, these terms are zero. When the 
evaluation point is at the position of the charge, the sum of the terms proportional to a delta function vanishes. From 
this expression the force is calculated according to the Lorentz force formula, specifically, = qPf^^u" . Here, the 
only non-zero component of the Maxwell field strength tensor is Frt = ^t,r, and the only non-zero component of the 
force is therefore the radial component. 



B. Regularization procedures 



1. Extended particle model: radial dumbbell 



Let us now assume for simplicity that the charges qi and 52 are separated only radially (there is no need for the more 
complicated splitting we did in the scalar case, because in the electric case we still need to perform mass regularization 
even for radial splitting). As before, we sum the forces at a common point at tq. After parallel transporting the forces 
radially to the common point tq, in the same way it was done above for the scalar case, we find that 
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(36) 



We do not regularize this bare force in the limit e ^ with the help of the exact solution because of the following. 
For an electric dumbbell in arbitrary acceleration in flat spacetime, the divergent piece of the self force is well known 
[p9[ , and is given by 



^div 
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(37) 



Here, a is the three acceleration, d is a unit three-vector in the direction of the dumbbell axis, and d is the length of 
the dumbbell axis. Note a factor of 2 between this expression and Eq. (57) of Ref. which is due to the fact that 
according to Eq. ( |l3|) the total force is twice the sum of the two forces. One would expect a similar expression to 



14 



hold also in curved spacetime. In our problem, the dumbbell axis is aligned in the radial direction, such that instead 
of a + ( a • d) 0? we would have 2or. (The acceleration is only radial.) We find 
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Note, that the length of the dumbbell axis is given by the invariant distance between ri and r2. As in the scalar case, 
we perform mass renormalization by subtracting this divergent piece of the force from the total force given by Eq. 
(^). Figure 1^ displays the renormalized force as a function of e. We find that the renormalized force approaches the 



correct finite value of f°^^^^ = e^M/(r^-y/l — 2M/r) |29y26t] like e, as indeed we expect 
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FIG. 7. Top panel: The renormalized force as a function of e. Bottom panel: |/^ 
is located at r = 2.1M. 
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2. Mode-sum regularization 



As in the scalar case, we can also construe the charge as point-like, and find from Eq. ( pq ) the total bare radial self 
force to be given by 
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For calculation of the bare force, Eq. ( p9| ) can be re-written as 
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which simphfies the calculation. However, Eq. (^0|) mixes the contributions of the different modes. Although the 
regularization procedure works also with this mixing, we shall consider below the regularization procedure with the 
force as given by Eq. (^9[). We first check the behavior of the modes '^^"/^(''o) as I oo. Figure ^ shows that 
indeed '^^"/^(^'o) approaches a constant, and that the difference between two consecutive modes scales like Z""^ for 
large values of in a similar way to the behavior of the modes for the scalar case. Consequently, also for this case, 
we infer that Ur ^ and that Cr — 0- We emphasize that for the case of an electric charge these parameters have not 
been calculated analytically, whereas in the scalar case they have. 



-0.95 




10° 10' 10^ 10^ 10* 

/ 

FIG. 8. Behavior of the I modes of the bare force for large values of /. Top panel (A): The behavior of ''^'^°F^(ro) as function 
of /. Bottom panel (B): | i^^(ro) - ''''■■°F^"^(ro)i as a function of /. The charge is located at r = 2.1M. 

In this case we also don't have prior knowledge about the value of the parameter dr. This is in general a serious 
problem, because without a knowledge of dr, the final result for the self force is not unambiguous. However, in this 
case we do have the final result from independent approaches, such that we can, in fact, predict the value of dr. It 
remains work yet to be done to compute the values of Or, br, Cr, and dr analytically for this case. 

As in the scalar case, we approximate the bare and the renormalized forces by the sum over a finite number of modes, 
and denote them by ^are^^ g^^^ ''"'^F^, respectively. We again define '^^'^F^ as in the scalar case, by subtracting 
from each mode of the series. Figure 9 shows the renormalized force ''''"F/ and its difference from Z™^'^* as functions of 
the mode number I. We find that ''''"F/ — f^^^'^^ approaches zero like for large values of I, such that we recover the 
results of Refs. i.e., we find that the self force is a repelling force, which is given by fr = q"^ M / {r^ ^yT^^2M/r) . 

The asymptotic agreement of and fr^'^^^ imply that also for this case dr = 0. This prediction can be tested 

analytically. 
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FIG. 9. Top panel (A): The renormalized force ^"''F' as a function of I. Bottom panel (B): \F['"^ - | as a function of /. 

The charge is located at r = 2.1M. The regularization procedure is performed with I' = 2.8 x 10^. 

IV. SUMMARY 

We presented a direct calculation of the self forces acting on two types of static charges in Schwarzschild spacetime: 
a scalar charge and an electric charge. In both cases the boundary conditions were chosen such that the scalar field and 
the potential, correspondingly, would be regular both at infinity and at the black hole's event horizon. Our method 
is based on decomposition of the field and the force into modes. Each mode satisfies an ordinary differential equation 
which we solve exactly in terms of Legendre Functions (in the scalar case) or derivatives of the Legendre functions 
(in the electric case). We find the total bare forces by summing over all modes numerically. This total force typically 
diverges. We then regularize the divergent self force with two independent procedures: First, we model the point-like 
particle to be spatially extended, and then consider a sequence of such particles, letting the spatial extention decrease. 
The divergent piece of the force is removed by a mass renormalization procedure (i.e., it is used to redefine the mass 
of the particle), and the remaining force approaches the self force in the limit of vanishing extension. Second, we use 
Ori's mode-sum regularization prescription, and remove the divergent piece of the force by studying the behavior of 
the bare force at large values of the mode number, and subtracting the value of the bare force at the limit of infinite 
mode numbers from all modes. Both regularization procedures recovered the well known results for static charges in 
the spacetime of a Schwarzschild black hole: a zero self force in the scalar case, and a repelling radial self force in the 
electric case. 

When one compares the relative effectiveness of the two regularization procedures, one finds that their effectivenesses 
are comparable. Specifically, for comparable values of I up to which we sum the series, we find that for both 
regularization schemes we obtain similar deviations of the computed regularized forces from the exact solutions, with 
roughly the same computation time. 

Evidently, more work is of need for both regularization prescriptions. In particular, it is not presently understood 
how to apply MSRP for more complicated cases, e.g., it is not presently clear whether there are cases with non- 
vanishing parameters c^, and whether the formalism can be extended to handle such cases (the radial component Cr 
was shown to be zero only for a scalar charge, although for all orbits in Schwarzschild). Also, it is not clear when 
non-zero functions should be expected. A generalization of MSRP to include also the gravitational case is also 
needed, a case for which the inherent gauge problem should be solved. We are currently using MSRP to study more 
complicated cases, in particular the self force acting on a scalar charge in circular orbit around a Schwarzschild black 
hole [|0|. 
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APPENDIX A: MODE-SUM REGULARIZATION PRESCRIPTION (MSRP) 



In this Appendix we describe very briefly the main ideas behind Ori's method for regularizing the mode-sum 



(MSRP) [g4 25| for a scalar charge in Schwarzschild. 

We emphasize that the work on this method is still in progress. However, for the case of a static scalar charge in 
Schwarzschild, the regularization scheme has been developed in full. 

As was pointed out by Quinn and Wald the physical self- force is the sum of two parts: (i) A local, Abraham- 
Lorentz-Dirac type term, and (ii) a "tail" term ^'^''F^, associated with the tail part of the Green's function. The local 
term is trivial to calculate (and it anyway vanishes in the static case considered in this paper). We shall therefore 
consider here the tail term only. This term may be expressed as 

= hm ,F^, (Al) 

where denotes the contribution to the force (evaluated at t = 0) from the part t < e of the particle's world line. 
Decomposing this expression into £-modes, one finds 



tail 



= }iT E ^//^ = }iZ E(''"/m - ^^f,)- (A2) 

Here, e/^, Seff^, and ^^"^/^ denote the force from the ^-multipole of the field sourccd by the interval r < e, the interval 
T > e, and the entire world line, respectively. The force ^^'"^ff^ may be identified with the sum over m and uj of the 
contributions from all stationary Teukolsky modes £, m, lu for a given £ (recall that in calculating a stationary field's 
mode £, m, uj one takes the source term to be the entire world line). Since we are using the retarded Green's function, 
the part r > does not contribute. However, the interval from e to 0+ does contribute. Essentially, it is this part 
which is responsible to the instantaneous, divergent, piece of the Green's function, which should be removed from the 
expression for ''^'i P 



A clarification is required here concerning the meaning of the last equality in Eq. (A2): Let tq denote the value of 
r at the evaluation point. Then, ^/^ is well-defined at r = rg. The situation with ^^^°ff^ and S^ff^ is more involved, 
however. Each of these quantities has a well-defined value at the limit r ^ r,^, and a well-defined value at the limit 
r Tq . Generically, for the r-component (and in some cases also for other components) these two one-sided values 



are not the same. Equation (A_2) should thus be viewed as an equation for either the limit r ^ oi the relevant 
quantities (i.e., ^^^°ff, and S^fJJ^or the limit r ^ of these quantities. Obviously, this equation is also valid for the 
averaged force, i.e., the average of these two one-sided values. In what follows we shall always consider the averaged 
force. Of course, the final result of the calculation, ^'^''F^ (which has a well-defined value at the evaluation point), is 
the same regardless of whether it is derived from its one-sided limit r — > r(^, or from r ^ r^, or from their average. 

Next, we seek an e-independent function /i^, such that the series J2ei^^''°ffi ~ ^ft) converges. Once such a function 
is found, then Eq. (|A2D becomes 



taii^^ = ^(^-7^ - h^) - £m Y^iSj'^ - hi). (A3) 



In principle, can be found by investigating the asymptotic behavior of ^'"'^ ff^ as £ ~> oo. It is also possible, however, 
to derive from the large-^ asymptotic behavior of Seff^ (the latter and ff^ must have the same large-^ asymptotic 
behavior, because their difference yields a convergent sum over £). In addition to ft.^, the investigation of S^f^ should 
also yield the parameter = lim^^Q- J2ii^tffi — h^), required for the calculation of ''^''F^ in Eq. (A3). 



Since we only need the asymptotic behavior of S^fj^ for arbitrarily small |e|, it is possible to analyze it using local 
analytic methods. In particular, we can apply a perturbation analysis to the £-TaoAe field equation (in the time 
domain). That is, we express the ^-mode effective potential V^{r) as a small perturbation 5V^{r) over the value 
of V'{r) at the evaluation point, Vq = V^[r = tq). Expressing G^[x^^ ,x'^{t)\, the £-mode Green's function, as a 
function of r and z = t£, the perturbation analysis provides an expression for G^lx'^ , x'^{t)] as a power series in r 
(with z-dependent coefficients). Only terms up to order are required for the calculation of the self force (recall 
that eventually we take the limit e ^ 0), and the perturbation analysis yields explicit expressions for the required 
three expansion coefficients of (as functions of z). Constructing Seff^ from G^ (essentially by integrating the 
latter's gradient from e to r = 0), it can be shown that the large-^ asymptotic behavior of S^ff^ takes the form 
^fiffi ^ "^^M^ + + c^£r^ + 0(^^^), in which the parameters 0^,6^,0^ are independent of £ and e (though they 
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depend on the orbit and evaluation point). (It can also be shown that there is no logarithmic divergence of ft.^.) The 
regularization function thus takes the form hj^ = a^£ + fo^ + Cfj.i~^, and the tail part of the self force is given by 

= E ("'"/m - - ^ " - d^. (A4) 

i 

In the case of a static scalar particle in Schwarzschild, one can show that = c,^ = = [p4| . (a^ and c^j are 
likely to vanish for all orbits in Schwarzschild, but so far this has been shown explicitly for the radial component 
only.) The self- force for a static particle then takes the simple form 



tail 



e 

Namely, in this simple case the regularization procedure is reduced to subtracting barcj£-^oo^ ^j^g large-£ limit of the i 
multipole of the bare force, from each multipole £ (note that since = 0, 6^ = barcj£^oo-j p^j. ^j^g particular case of 
a static scalar charge in Schwarzschild, Ori j2j] also obtained analytically the value of this largc-i? limit of the force: 
b^^-[qy{2r^)]{l-M/r)/il-2M/r)S;. 

This regularization prescription takes a trivial form in the cases of static scalar or electric charges in Minkowski 
spacetime. In these cases it is easy to verify that all the £ modes of the bare force are equal (i.e., independent of £), 
specifically — / (2r^) = const (this can be obtained easily directly from a decomposition of the field). When 

one sums over all modes, the bare force of course diverges. However, subtracting this constant term from each mode 
yields a new series, where all modes are zero, such that the total force vanishes, which is the well-known result in 
Minkowski spacetime. We note that MSRP turns out to be effective also for the cases of scalar or electric charges in 
circular orbits in Minkowski spacetime . 

We emphasize that whereas the parameters a^, 6^, and can be found from the behavior of ff^ at large values 
of the mode number £, the parameter can only be calculated according to its definition. In the simple case of 
a scalar charge in circular orbit around Schwarzschild, which includes as a special case a static scalar charge, this 
calculation is not hard to do, and the exact value of was found (for this case = 0). However, it might be the 
case that for more complicated cases d^ is more difhcult to find. Then, one can still regularize the force, but the 
result would not be free from the ambiguity which results from the ignorance of the exact value of d^. 

MSRP involves integration over the entire world line of the orbiting object. In that respect, it is especially suitable 
for periodic, or near periodic, orbits. For a-periodic orbits, such as the final plunge of the object into the black hole, 
one can perhaps use a different approach, where one integrates only over the past world line, excluding the position 
of the particle itself, and thus avoids the singular contribution. Of course, the closer to the particle one integrates, 
the more modes one would need to sum over in order to obtain convergence. In fact, the number of modes is inversely 
proportional to the proper-time difference from the event up to which one integrates and the position of the particle. 
It has been recently shown by Wiseman that in the far-field limit (i.e., to leading order in the ratio of the black 
hole mass to the radius of the orbit), the contribution of the near neighborhood of the past world line is negligible, 
such that one may need to sum only over a relatively small number of modes. In stronger fields, this approach can 
perhaps be combined with a normal neighborhood expansion 0,0 to obtain the self force. 



APPENDIX B: NUMERICAL EVALUATION OF THE LEGENDRE FUNCTIONS 

All the series we need to evaluate very accurately involve the Legendre functions of the first and second kinds, and 
their derivatives. The degree of the functions is very high. For example, in the numerical summations reported here 
we evaluated the series up to Z = 4.5 x 10*. This requires us to use a very accurate algorithm for the calculation of 
the Legendre functions. In fact, we find that the Legendre functions of the first kind can be very accurately computed 
by using the recursion relations, such as 

Pi+i{x) = [{21 + l)xPi{x) - IPi^iix)] (Bl) 

and 
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Although similar relations hold also for the Legendre function of the second kind, they are not practical for the 
following reason. The functions Qi{x) approach zero very fast for fixed x > 1 when the degree gets very large. The 
subtraction which is inherent to the recursive expression becomes numerically inaccurate very rapidly. The functions 
Qi can also be considered as the sum of two series, one being a polynomial, and the other being a polynomial multiplied 
by a common logarithmic factor. Each of the polynomials satisfy the same recursive formula as the Legendre functions, 
but with different initial terms for I = and 1 = 1. Each of the two series grow very fast with /, but their difference 
gets very small. Therefore, this method would also not be very accurate numerically. A way to avoid these difficulties 
is to use the integral representation of the functions Qi{x). This is given by 

The integrand does not have any pathologies in the entire interval of integration, and also the boundaries are regular. 
We perform this integral using Romberg integration, which proves to be very efficient and accurate p3| . The derivatives 
of the functions Qi{x) can still be computed by the relation given above for Pl{x). Another improvement on the 
numerical evaluation of both Pi{x), Qi{x) and their derivatives is the following. In all the expressions we have, we 
need only compute the product of two Legendre functions or their derivatives, one factor involving Pi (or its derivative), 
and the other involving Qi (or its derivative). Because we are not interested in the value of the Legendre functions 
themselves, but only in such products, we can disregard the factor of 2~^'+^^ in the integral representation of the 
Qi{x). This would mean that each of the functions we compute is too large by a factor of 2'+^. If we then compute, 
instead of the functions Pi{x), a new function, which is smaller than Pi{x) by the same factor, the product of the 
two new functions would be unchanged. This can be done also for the derivatives of the Legendre functions. It is 
advantageous to do this, because for a given floating point arithmetic this procedure increases the maximal value of 
I for which accurate computations can be performed by an order of magnitude. 
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